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Abstract. The Hudson Bay Lowlands (HBL) is the second 
largest boreal wetland ecosystem in the world and an impor- 
tant natural source of global atmospheric methane. We quan- 
tify the HBL methane emissions by using the GEOS-Chem 
chemical transport model to simulate aircraft measurements 
over the HBL from the ARCTAS and pre-HIPPO campaigns 
in May-July 2008, together with continuous 2004-2008 sur- 
face observations at Eraserdale (southern edge of HBL) and 
Alert (Arctic background). The difference in methane con- 
centrations between Eraserdale and Alert is shown to be a 
good indicator of HBL emissions, and implies a sharp sea- 
sonal onset of emissions in late May (consistent with the air- 
craft data), a peak in July-August, and a seasonal shut-off in 
September. The model, in which seasonal variation of emis- 
sion is mainly driven by surface temperature, reproduces well 
the observations in summer but its seasonal shoulders are too 
broad. We suggest that this reflects the suppression of emis- 
sions by snow cover and greatly improve the model simu- 
lation by accounting for this effect. Our resulting best esti- 
mate for HBL methane emissions is 2.3 Tg a“^ several-fold 
higher than previous estimates (Roulet et al., 1994; Worthy 
et al., 2000). 
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1 Introduction 

Methane is the second most important anthropogenic green- 
house gas after carbon dioxide (IPCC, 2007). Methane 
concentrations have increased from 700 ppbv in the pre- 
industrial atmosphere to 1700 ppbv by the early 1990s 
(Etheridge et al., 1998). This increase is presumably driven 
by direct emissions from industry and agriculture (IPCC, 
2007), but could also reflect changes in the chemical sink (re- 
action with the OH radical) and the effects of climate change 
on natural emissions (Worthy et al., 2000). Wetlands are the 
largest natural source of methane and are highly sensitive to 
changes in climate (Kaplan et al., 2006), especially in the bo- 
real zone (Zhuang et al., 2006; Sitch et al., 2007). Here we 
use aircraft observations over the Hudson Bay Lowlands in 
northern Ontario as well as surface observations at Eraserdale 
and Alert to better quantify this boreal wetland source. 

Methane is produced in boreal wetlands by heterotrophic 
respiration of partially decomposed organic material under 
anoxic conditions. As it rises to the surface, it may en- 
counter oxic conditions resulting in oxidation. The remain- 
ing methane escapes to the atmosphere (Walter et al., 2001). 
Production of methane primarily depends on soil tempera- 
ture, and the distribution of anoxic and oxic zones depends 
primarily on the water table level (Pelletier et al., 2007; 
Moore et al., 1994). Compounding this complexity in bo- 
real wetland emissions is the seasonal thaw that propagates 
from the surface to depth (Dunn et al., 2009; Wania et al., 
2009). 
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The Hudson Bay Lowlands (HBL) is an ecologically 
significant and well- studied boreal wetlands region (e.g. 
Glooschenko et ah, 1994). It is (after the West Siberian wet- 
lands) the second largest semi-continuous wetland region in 
the world, covering an area of 320 000 km^ or about 10% of 
the total area covered by boreal wetlands (Wang et ah, 2008; 
Glooschenko et ah, 1994). The ABLE-3B/NOWES airborne 
and ground campaign conducted in the summer of 1990 esti- 
mated an annual methane emission of 0.5 zb 0.3 Tg a“^ from 
the HBL (Roulet et al., 1994). Worthy et al. (2000), using 
inverse methods to interpret observations from the Alert and 
Eraserdale Canadian sites, obtained a similar estimate of 0.2- 
0.5Tga“^ These estimates seem low compared to global 
inversion studies that infer total boreal wetland emissions of 
27-38 Tga“^ (Hein et al., 1997; Wang et al., 2004; Bous- 
quet et al., 2006), although Chen and Prinn (2006) obtained 
an estimate of 7Tga“^. Typically, the HBL is assumed to 
contribute 10% to global boreal wetland emissions. 

The above inconsistencies point to the need for a bet- 
ter understanding of methane emissions from the HBL as a 
window into the global boreal wetland source. We exploit 
here methane concentration measurements from the ARC- 
TAS and Pre-HIPPO aircraft campaigns in May-July 2008 
(Pig. 1), together with long-term surface data from Envi- 
ronment Canada at Eraserdale (81.6° W, 49.9° N) and Alert 
(62.5° W, 82.5° N). We interpret these data using a global 
bottom-up scheme for wetland emissions implemented in the 
GEOS-Chem chemical transport model (CTM). The ARC- 
TAS airborne campaign based in Cold Lake, Alberta (Jacob 
et al., 2010) conducted three flights over the HBL in early 
July, while the Pre-HIPPO campaign based in Boulder, Col- 
orado (Pan et al., 2010) conducted 2 flights in the region in 
May-June. 


2 Model description 

We use the GEOS-Chem CTM originally described by Bey 
et al. (2001) in a methane simulation for 2004-2008 to in- 
terpret the aircraft and surface observations. GEOS-Chem 
is driven by GEOS -5 analyzed meteorological data from the 
NASA Global Modeling and Assimilation Office (GMAO). 
The GEOS-5 data have 1/2° x 2/3° (lat. Ion) horizontal reso- 
lution with 72 vertical levels and 6-h temporal resolution (3- 
h for surface variables and mixing depths). The horizontal 
resolution is degraded here to 2° x 2.5° for input to GEOS- 
Chem. 

The methane simulation in GEOS-Chem was originally 
described by Wang et al. (2004), and subsequently improved 
and updated by Drevet et al. (2011). Major sources in- 
clude anthropogenic emissions from EDGAR 4.0 (Euro- 
pean Commission, 2009), natural emissions from wetlands 
as described below, and GPED2 biomass burning emissions 
(Giglio et al., 2006). Chemical loss of methane is com- 
puted using a global 3-D archive of monthly average OH 
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Fig. 1. ARCTAS and Pre-HIPPO flight tracks over the Hudson Bay 
Lowlands (HBL) below 4 km, superimposed on a map of GEOS- 
Chem methane emissions for July 2008. The locations of Eraserdale 
and of the two ABLE-3B/NOWES study regions are also shown. 
The black rectangle encompasses the HBL region as defined in the 
present study. 


concentrations from a GEOS-Chem simulation of tropo- 
spheric chemistry. The mean tropospheric OH concentra- 
tion is 10.8 X 10^ molecules cm“^, typical of current models 
(a global model inter-comparison by Shindell et al. (2006) 
gives a mean value of 11.1 ± 1.7 x 10^ molecules cm“^). 
The corresponding tropospheric lifetime of methylchloro- 
form against oxidation by OH is 5.3 yr, within the range con- 
strained by observations (Prinn et al., 2005). Additional mi- 
nor sinks for methane in the model include stratospheric oxi- 
dation prescribed as a constant decay (Wang et al., 2004) and 
soil absorption (from EDGAR 4.0). 

The wetlands emission scheme in GEOS-Chem is based 
on Kaplan (2002). The scheme applies methane emission 
factors to heterotrophic respiration rates in tropical and bo- 
real wetlands, following algorithms described by Kaplan et 
al. (2002), Sitch et al. (2003) and Bergamaschi et al. (2007). 
The emission flux E (molecules CH 4 m“^s“^) for each 
model grid square is given by: 


-Eq 

E = 8WFfiAe'^-^o 


E 


Q 


Ei 


( 1 ) 


where C\ and C 2 (molCm“^) are soil and litter carbon 
pools, respectively, specified on a 2° x 2.5° grid by the Lund- 
Potsdam-Jena dynamic global vegetation model (Sitch et al., 
2003). These pools have fixed residence times ri = 32 yr 
and Ti = 2.8 yr from Sitch et al. (2003). The Arrhenius fac- 
tor with A = 1.0^+^ Eo = 309 K, and Tq = 221 K (Lloyd 
and Taylor, 1994) specifies the temperature dependence of 
respiration. The soil temperature T is taken here as the 
GEOS-5 skin temperature. The methane emission factor is 
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P = 3 X 10“^ mol CH4 per mol C respired (Christensen et 
ah, 1996). An additional scaling factor F is used to match 
observed ecosystem fluxes of methane, separately for tropi- 
cal (T) and boreal (B) wetlands: 

F = a • Fj -h ( 1 — O') • F^ (2) 

where a = min(exp[(7' — ri)/8] , l), f (K) is the mean 
monthly soil temperature taken here as the GEOS-5 skin tem- 
perature, and T\ = 303.15 K. Scaling factors Fj = 0.14 and 
Fb = 0.005 were derived by J. Kaplan and J. Drevet (per- 
sonal communication, 2009) to match published emission 
estimates for the Amazon (Melack et al., 2004) and boreal 
wetlands (Hein et al., 1997; Wang et al., 2004). Finally, 
W in Eq. (1) represents the maximum potential fraction of 
wetland coverage for the grid square, as obtained from vari- 
ous databases described by Kaplan et al. (2002) and Berga- 
maschi et al. (2007). Whether a wetland is actually present 
(8 = 1) or not (8 = 0) over that fraction W at any given time 
depends on the level of the water table, which we diagnose 
using GEOS-5 soil moisture as a proxy following Bergam- 
aschi et al. (2007). Wetlands are present if the soil moisture 
exceeds a specifled threshold of 0. 1 for the ratio between the 
soil water content and the porosity of the soil (J. O. K. Ka- 
plan, personal communication, 2009) 

Annual emissions for the HBL (geographically deflned as 
50° N-60° N, 75° W-96° W) computed in the above manner 
in GEOS-Chem average 2.9Tga“^ for 2004-2008, with the 
spatial distribution shown in Fig. 1. This estimate is much 
larger than the previous HBL emission estimates of Roulet et 
al. (1994) and Worthy et al. (2000). 

3 Constraints on HBL methane emissions 

We examine the consistency of the GEOS-Chem methane 
emissions with the aircraft data from Pre-HIPPO and ARC- 
TAS, and the surface data from Fraserdale and Alert. The 
DACOM tunable diode laser instrument used in ARCTAS 
(Sachse et al., 1987) has an estimated accuracy/precision 
of 1%/0.1%. The quantum cascade laser instrument used 
in Pre-HIPPO has an estimated accuracy and precision of 
0.25%. The surface measurements at Fraserdale are ob- 
tained by gas chromatography on samples collected from a 
40-m high tower and have an accuracy/ precision of 1 %/0.2% 
(Worthy et al., 2003). Similar speciflcations apply to the sur- 
face measurements at Alert. 

Figure 2 shows the ensemble of aircraft vertical profiles 
over the HBL from 12 May (Pre-HIPPO) to 5 July (ARC- 
TAS). We excluded stratospheric air as diagnosed by a molar 
O3/CO ratio exceeding 1.25 (Hudman et al., 2006) and Are 
plumes as diagnosed by CO exceeding 200ppbv. The lat- 
ter Alter effectively removed biomass burning influence from 
the data set as inferred from correlation between methane and 
CO. The HBL methane enhancements in Fig. 2 can thus be 
reliably attributed to wetland emissions. 


The ARCTAS observations on 4-5 July show strong 
boundary layer enhancements over the HBL. The Pre-HIPPO 
flight on 12 May shows no boundary layer enhancement 
while that on 23 June shows a moderate enhancement. Ob- 
servers on the Pre-HIPPO aircraft reported snow cover over 
the HBL on 12 May but not on 23 June. For comparison 
with the aircraft we sample the model at the time and loca- 
tion of the flights. We see in Fig. 2 that the model provides 
a good simulation of the boundary layer structure for the dif- 
ferent flights, the enhancement observed in ARCTAS, and 
the sharp springtime transition from May to July. However, 
model overestimation is evident for the 23 June profile. 

To further investigate the magnitude and seasonal on- 
set of HBL emissions we used 2004-2008 surface data at 
Fraserdale and Alert collected by Environment Canada, with 
Alert serving as an Arctic background site against which the 
HBL influence at the Fraserdale downwind site can be refer- 
enced (Worthy et al., 1998). For Fraserdale we sample the 
daily data averaged over the 1700-1900 local time window, 
when the surface measurements are most representative of a 
relatively deep mixed layer (Worthy et al., 1998), and further 
select for surface winds from the northern quadrants, when 
direct influence from the HBL can be expected (Fig. 1). Se- 
lection for northerly winds retains ~50% of the original data. 
We sample Alert data for the same times at Fraserdale in or- 
der to facilitate analysis of the difference between the two 
sites as discussed later. 

Figure 3 shows the observed seasonal variations at 
Fraserdale and Alert for 2004-2008. The observations at 
Alert show a July minimum due to chemical loss in the 
Northern Hemisphere. The model minimum lags 4-6 weeks 
behind, an offset that can be attributed to background error 
in the seasonal variation of sources, transport, or OH concen- 
trations. The observations at Fraserdale follow the seasonal 
variation at Alert in winter- spring but deviate in late May 
toward an August maximum, ostensibly due to emissions 
from the HBL. The model shows the same seasonal devia- 
tion at Fraserdale relative to Alert but shifted 6 weeks early. 
A model sensitivity simulation with no HBL emissions (also 
shown in Fig. 3) confirms that the seasonal deviation between 
Fraserdale and Alert is due primarily to HBL emissions. 
The model shows multiple seasonal peaks at Fraserdale (late 
June, late August, early November) compared to a single ob- 
served peak, but this fine structure reflects fluctuations in the 
background rather than HBL emissions as discussed below. 

Although Fraserdale is at the southern end of the HBL 
and the data are collected at only 40 m altitude, they ap- 
pear to be reasonably representative of the HBL. As a test, 
we partitioned the HBL into northern and southern halves 
and examined the model sensitivity of Fraserdale to emis- 
sions from each half. We found no significant difference. 
We also found no significant difference in methane concen- 
trations when sampling the model at different altitudes over 
Fraserdale up to 500 m (five model layers). 
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Pre-HIPPO 


ARCTAS 


Fig. 2. Methane vertical profiles from Pre-HIPPO and ARCTAS over the HBL (May-July 2008). Observations (blue) are compared to 
GEOS-Chem (GC) model vertical profiles sampled along the flight tracks at the flight times. The standard simulation (red) and a sensitivity 
simulation with no HBL emissions (black) are presented. 



Fig. 3. Seasonal variation (2004-2008) of methane at Fraserdale 
and Alert. GEOS-Chem results are compared to observations. Also 
plotted is the model background concentration at Fraserdale as de- 
rived from a simulation with no HBL emissions. Data are daytime 
values smoothed with a 28 -day moving average and then averaged 
over 5 yr. For Fraserdale we use only data associated with winds 
from the northern quadrants. 

Figure 4 shows the seasonal variation of the difference in 
concentrations between Fraserdale and Alert (ACH 4 ), illus- 
trating more precisely the methane flux signature from the 
HBL. Here we assume that Alert provides a reasonable mea- 
sure of background concentrations at Fraserdale; this is sup- 
ported in the model by the comparison in Fig. 3 of the model 
simulation at Alert (thin blue line) and at Fraserdale in the ab- 
sence of HBL emissions (thick black line). We And that the 
multi-peak structure of model concentrations at Fraserdale 
in June-November (Fig. 3) is reduced when corrected for the 



Fig. 4. Mean seasonal differences in CH 4 concentration between 
Fraserdale and Alert (ACH 4 ) for 2004-2008 (data in Fig. 3). Ob- 
servations (blue) are compared to the standard GEOS-Chem sim- 
ulation, a sensitivity simulation restricting emissions to snow-free 
ground, and a sensitivity simulation with no HBL emissions. 

Alert background (Fig. 4). More importantly, we And that 
any residual multi-peaks in modeled ACH 4 are associated 
more with changes in the model background at Fraserdale 
relative to that at Alert than in HBL model emissions. Tem- 
poral fluctuations during those two months in the model may 
also reflect the greater variability in surface temperatures 
(used in the model to compute methane emission) than in ac- 
tual soil temperatures. Heat transfer in the soil column would 
be expected to dampen temporal variability in soil tempera- 
tures. 

The observed onset of HBL methane emissions in late 
May, as seen in Fig. 4, is consistent with the aircraft 
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observations discussed above and with previous field stud- 
ies in nearby James Bay peatlands that suggest an onset of 
emissions in mid-May (Pelletier et al., 2007). By contrast, 
we see from Fig. 4 that HBL emissions in the model begin 
in early April. In addition, the observations indicate a sea- 
sonal shutdown of HBL emissions in September whereas in 
the model these emissions persist into October. The early on- 
set of model emissions was not apparent in comparison with 
the Pre-HIPPO profile on 12 May in Fig. 2, but that is be- 
cause of delayed spring warming in 2008 and because that 
flight profile sampled the northern edge of the HBL (Fig. 1). 

The premature onset of HBL methane emissions in the 
model likely refiects the use of skin temperature as proxy 
for soil temperature. Seasonal increases in soil temperature 
at depth lag behind the land surface during the spring thaw. 
We attempted to impose in the model a time lag for soil heat- 
ing by using the standard heat transport parameterization of 
Campbell and Norman (1998) with thermal diffusivities from 
Sitch et al. (2003), but the resulting delay in the onset of 
emissions was insufficient. Instead we identified persistent 
snow cover in the GEOS-5 data well past the model onset in 
model emissions. Snow cover would insulate the underlying 
soil from warming, inhibiting methanogenesis in spring, and 
would also trap methane in the autumn (Friborg et al., 1997). 
Consequently we modified the model to restrict emissions to 
snow-free regions. Figure 4 shows that this mostly corrects 
the model biases in the spring and autumn, although there is 
still a small lag of 1-2 weeks in the spring. This additional 
delay might refiect a period of time required for the underly- 
ing peatlands to thaw before methanogenesis ensues (Dunn 
et al., 2009). Comparison with the aircraft profiles in Fig. 2 
is unaffected by the delayed onset in model emissions. The 
resulting annual reduction in model HBL emissions is 20% 
(2.3 Tg a-* vs. 2.9 Tg a-*). 

The model temporal variability of methane emissions in 
the snow-free season is driven largely by surface tempera- 
ture (Eq. 1), and this appears adequate to match the observed 
July-August maximum of HBL emissions (Fig. 4). Previous 
studies of boreal wetlands have pointed out the sensitivity of 
emissions to changes in the level of the water table (Moore et 
al., 1994; Pelletier et al., 2007). However, the fiat topography 
of the HBL results in poor drainage and maintains persistent 
wetland coverage throughout the summer. 

Figure 5 shows the interannual variation of HBL model 
emissions for 2004-2008 as driven by temperature and snow 
cover. The seasonal onset of emission can vary by a month 
from year to year. There is much less year-to-year variability 
in the fall shutdown of emissions. The mean annual emission 
for the 5 yr is 2.3 zb 0.3 Tg a“^ . 

4 Comparison to ABLE-3B/NOWES estimates 

The ABLE-3B/NOWES surface and aircraft field study in 
July 1990 previously reported an annual emission estimate 



Fig. 5. Seasonal variation of HBL methane emissions simulated 
by the model for 2004-2008. Values are integrated spatially over 
the HBL domain (50° N-60° N, 75° W-96° W) and smoothed tem- 
porally with a 4-week moving average. Also tabulated are annual 
emission estimates for individual years. 

of 0.5 di0.2Tga“^ for the HBL (Harris et al., 1994; Roulet 
et al., 1994). This is considerably less than our best esti- 
mate of 2.3 Tg a“^ and would be inconsistent with the Pre- 
HIPPO and ARCTAS data of Fig. 2 as well as the Fraserdale 
ACH4 data of Fig. 4. The ABLE-3B/NOWES estimate was 
obtained by extrapolation of direct fiux measurements at sur- 
face sites, using wetland coverage derived from satellite and 
aerial imagery. The surface sites and supporting aircraft eddy 
correlation fiux measurements were located in two small 
study areas at the southern and northern edges of the HBL 
(Fig. 1). Roulet et al. (1994) reported mean June-October 
emission estimates for the southern and northern study ar- 
eas of 3.4gm“^a“^ and 6.3gm“^a“^ respectively, from 
their surface measurements. Aircraft measurements over 
these same regions in July yielded consistent mean fiuxes 
of 5 zb 3 g m~^ a“^ and 4 zb 6 g m“^ a“^ respectively. These 
values agree with our fiux estimates of ~5gm“^a“^ for 
both ABLE-3B/NOWES regions in July (Fig. 1). Roulet et 
al. (1994) went on to infer annual mean emissions by tak- 
ing their June-October measurements to be representative of 
snow-free conditions and assuming zero emissions for snow- 
covered ground. This is consistent with our findings. 

The large difference between our estimate of HBL emis- 
sion estimate and that of Roulet et al. (1994) thus lies in the 
spatial extrapolation to the scale of the HBL. We see from 
Fig. 1 that our emissions are much higher in the mid- section 
of the HBL than in the ABLE-3B/NOWES study regions. 
The boundary layer methane enhancements observed from 
the ABLE-3B aircraft ("^30ppbv) were indeed much lower 
than the mean value of 60ppbv observed on the ARCTAS 
flights (Fig. 2). 
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5 Conclusions 

Aircraft observations over the Hudson Bay Lowlands (HBL) 
in May-July 2008 show a seasonal onset of methane emis- 
sions in June and 60ppbv enhancements in the boundary 
layer in July. Surface observations at Fraserdale (just south 
of the HBL), when referenced against a background Arctic 
site (Alert) to isolate the HBL contribution, indicate a sea- 
sonal onset of methane emission in late May, a peak emis- 
sion from mid-July to the end of August, and a sharp de- 
crease in September. The GEOS-Chem model including a 
standard methane emission scheme for boreal wetlands can 
successfully reproduce these observations except for a pre- 
mature springtime onset and a delayed fall shut-off. Seasonal 
variation of wetland emission in the model is mainly driven 
by surface temperature. We find that accounting in addition 
for suppression of emission by snow cover corrects the model 
biases in spring and fall. The variability in the model is still 
larger than observed and this could refiect dampening of soil 
temperature fluctuations relative to the surface. Our result- 
ing best estimate of HBL methane emissions is 2.3Tga“\ 
much higher than previous estimates for the region (Roulet 
et al., 1994; Worthy et al., 2000). We argue that this reflects 
gradients of methane emission within the HBL that were not 
previously accounted for. 
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